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Abstract 

We find exact solutions in five dimensional inhomogeneous matter dominated 
model with a varying cosmological constant. Adjusting arbitrary constants of in- 
tegration one can also achieve acceleration in our model. Aside from an initial 
singularity our spacetime is regular everywhere including the centre of the inhomo- 
geneous distribution. We also study the analogous homogeneous universe in (4+d) 
dimensions. Here an initially decelerating model is found to give late acceleration 
in conformity with the current observational demands. We also find that both 
anisotropy and number of dimensions have a role to play in determining the time 
of flip, in fact the flip is delayed in multidimensional models. Some astrophysical 
parameters like the age, luminosity distance etc are also calculated and the influence 
of extra dimensions is briefly discussed. Interestingly our model yields a larger age 
of the universe compared to many other quintessential models. 
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Recent observational evidences suggest that the present day universe has the critical 
energy density containing presumably 70% dark energy and about 30% dark matter, where 
the term dark indicates a sort of invisibility. While the sceptics will always question 
the wisdom to explain data based on something we can not see the avalanche of data 
emanating from type la supernovae measurements [I], CMB anisotropics [2], galactic 
rotation curves and surveys of galaxies, clusters and superclusters make the presence of 
dark matter and dark energy increasingly convincing. This situation inevitably forces us 
to pose the question : why are dark matter and dark energy so dark ? 

As a natural corollary to the question one is reminded of another 'dark' stuff in physics, 
extra dimensions. Like dark matter and dark energy to explain the current quintessential 
behaviour the existence of extra dimensions is absolutely necessary in any attempt to 
unify gravity with other forces of nature beyond the standard model of particle physics 
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and also to explain the vexed hierarchial problem of quantum mechanics in varied brane 
inspired models [3j. The recent spurt in activities in extra dimensions also stems from the 
Space-Time-Matter (STM) theory [4] proposed recently by Wesson and his collaborators 
[5]. As our space time is manifestly four dimensional in nature the extra dimensions 
should be 'hidden' ( or 'dark' ). 

In the back drop of the above discussions it is more than apparent that this common 
feature of 'invisible existence' of dark matter, dark energy and extra dimension points 
to the conjecture that there must be some deep underlying relationship and inter play 
among them [6]. 

Motivated by this consideration we have in the past worked out problems [5] where 
starting from a perfect fluid in higher dimensions obeying reasonable energy conditions 
accelerating model is made possible as a consequence of the existence of higher dimensions. 
Interestingly we do not have to invoke any adhoc quintessential type of fluid by hand to 
achieve the acceleration. 

The motivation for the present work is somewhat different. In an earlier communica- 
tion [7] we have analysed an inhomogeneous 5D spacetime with a constant A, and showed 
that unlike the finding of Tosa [8] the spacetime did show the desirable feature of dimen- 
sional reduction as the usual 3D space expands. Although in the light of observational 
evidences that our universe is currently accelerating, the idea of A is now an integral in- 
gredient in cosmological models its introduction, nevertheless, invites serious conceptual 
problems. The upper limit of A from observations is about 120 orders of magnitude below 
the value for the vacuum energy density predicted from quantum field theory. 

To circumvent this difficulty and a host others an essentially phenomenological ap- 
proach is generally taken where it is argued that due to the coupling of the dynamical 
degrees of freedom with matter fields of the universe, A relaxes to its current small value 
through the expansion of the Universe and creation of photons [9]. From this point of 
view A is small because the Universe is old. 

Another salient feature of our model is its inhomogeneity. While the higher dimen- 
sional generalisation of the FRW models has been adequately addressed in the literature 
scant attention has been paid so far to address the inhomogeneous situations and the 
issues coming out from it. Moreover the work of Mustapha etal [TD] and others indi- 
cate that there is no unquestionable observational evidences for spatial homogeneity. So 
investigations in inhomogeneous models is always a welcome step. 

Although the cosmological concordance ACDM models fit most of the observations 
well, as mentioned earlier, it is beset with theoretical difficulties. In fact whenever there 
is any ambiguity interpreting any observational results as in 1990s regarding the age of 
the universe [11] the inhomogeneous models are seriously considered. Similarly in the 
backdrop of the current accelerating phase of the universe the role of inhomogeneities 
[T2] like that of extra dimensions as a possible ingredient of the cause of the apparent 
acceleration is being increasingly debated [13] . In fact averaging out the spatial degrees of 
freedom ( i.e. inhomogeneities) leads to averaged, effective Einstein equations (somewhat 
similar to what we find in higher dimensional models also), which contain some additional 
terms relating to the inhomogeneous distributions generally termed back reaction in this 
context. It is argued that this so called back reaction can account for the current accel- 
eration of the universe [14J. Our scope in the present work is very limited. While we 
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have attempted to address the implications of extra dimensions in the context of current 
accelerating expansion in some detail the role inhomogeneities in our present model has 
been deferred to a future work in progress. 

The paper is organised as follows: in section 2 we find an exact solution for a 5D 
inhomogeneous distribution with a variable A. Interestingly we here get acceleration, 
although at the cost of sacrificing dimensional reduction. In section 3 we focus on a 
homogeneous (d+4)-dimensional spacetime to show that accelerating model is possible in 
this case also. In section 4 some important astrophysical parameters are discussed for our 
model to investigate how inclusion of extra dimensions influences the situation. Section 
5 compares the A varying case with our earlier A- constant situation. Here an initially 
decelerating universe ends up as an accelerating one. Here we find two interesting results- 
both anisotropy and number of dimensions have marked effects in determining the instant 
of flip. The paper ends with a brief discussion in section 6. 



2. The Field Equations and its integrals 

We consider a metric of (3 + d + l)-dimensional spacetime where both the ordinary 
3-space and the extra ( or hidden) space are inhomogeneous, isotropic and flat : 

ds 2 = dt 2 - a 2 [dr 2 + r 2 {d6 2 + r 2 sin 2 6d<p 2 )] - b 2 dy 2 (1) 

dy 2 = Yli=i dyf and we assume a single scale factor for the internal dimensions, here a(r,t) 
and b(r,t) are scale factors. Assuming that the matter content in this higher dimensional 
space is taken to be a perfect fluid and is augmented by the inclusion of a time vary- 
ing cosmological constant we can write down the Einstein's equations, which govern the 
evolution of the ordinary 3-space as well as the extra space as 
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where the (3+d+l)-dimensional stress-tensor, Ty, in comoving co-ordinate should be of 
the form 

T ° = p + A, T\ = If = T| = -p + A, T d d = — p,j + A (7) 

We have here assumed that 8ttG = 1 where G is the gravitational constant in the higher 
dimensional space and a dot and and a prime overhead denote differentiation w.r.t. time 
and radial co-ordinate. 

To make the field equations tractable we assume, at this stage, that the scale factors 
a = a(t) and b = b(r,t) and d = 1, i.e., we take a 5D spacetime where the y = constant 
hypersurface is homogeneous FRW-like. In what follows we shall see that the fact that 
inhomogeneity is being introduced through the extra space has far reaching implications 
in the cosmological evolution of our models. For simplicity as also the fact that the present 
universe is matter dominated, we consider here the dust case only (p = p^ = 0). 

Now the equation (2) yields 

b = af3(r) + a(t) (8) 

where (3{r) and a(t) are integration functions of r and t respectively. Using equations 
(3), (4) and (8) we get 

b = - C r 2 a + a{t) (9) 

where c is a pure constant. However we subsequently see that c is not exactly arbitrary, 
being a measure of the curvature of the 4D space. Moreover c = implies at once 
a flat and homogeneous spacetime. In line with the existing literature one assumes the 
dynamical behaviour of A as A ~ (^) 2 , A ~ - or A ~ p. Although not strictly independent 
the first one was proposed from dimensional consideration by Lima and Waga [15] and 
subsequently adopted by several workers [16], [17] whereas the second type by Overduin 
and Cooperstock [18] while Vishwakarma [9] favoured the third alternative. 

3 

In our work we take A = kH 2 ,(H = ^) which through equation (5) yields a = t^~ k ) 
and A = ^ 6 ^ 2 t~ 2 . Now using equation (3) we get, 

6 i (3 + k) 9 a 

& + 7? — t\ t ~ 3 77 — r&t a + 4ct (6 " fc) = 10 
(o — k) (o — ky 

which through a long but somewhat straight forward calculation yields 

3 (3+fc) 2c(6 — k) 2 (9-2*0 , , 

a = ait «-«> + a 2t «-» - {12 \ ){3 '_ k) t «-« (ID 
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di, a 2 are arbitrary constants. So clubbing everything together we finally get 
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One can now calculate the 4-space curvature (i?*( 4 )) of the t-constant slice for the line 
element (l)(d = 1) through the expression, R\ = R*^ + 9 + 9 2 - 2cu 2 + [15] . where 9 
is the expansion scalar and last two terms refer to vorticity and acceleration. After some 
algebra we get the expression 

iT (4) = — (16) 

( fll - cr 2 )tl^) + a 2 r T^T - ffggg) 

So as commented earlier, the arbitrary constant c comes out to be a measure of the 
curvature of the 4D space, c is also a measure of inhomogeneity parameter because 
c = makes our spacetime homogeneous. Relevant to mention that unlike the analogous 
homogeneous case here the curvature also depends on spatial co-ordinate r. Moreover it 
blows off at the initial singularity t = 0. Otherwise it is regular everywhere including the 
origin of the distribution. 

To check if our space time contains any geometric singularity aside from the well known 
t — big bang epoch we calculate the Kretschmann scalar for our 5D inhomogeneous line 
element as 

+ + (17) 
a 2 b 2 a 4 \ab a 2 br I 

At £ — >■ 0, the invariant diverges but it is regular at r = since — is regular there. So 
unlike similar cases in many inhomogeneous distributions there is no spatial singularity in 
our cosmological model including the centre of distribution. Moreover it has not escaped 
our notice that while the 3D scale factor a starts from t = the 5D scale b is infinite 
there. To ensure that both the scales start simultaneously at the initial singularity we set 
the arbitrary constant a 2 to be zero henceforth. 



The inhomogeneous metric we presented here is very general in nature in the sense 
that many well known solutions in this field are recoverable as special cases of our metric. 

(i) When k = the A term vanishes and we get back our earlier solution [20] . 

(ii) If k = , c = ai = 0, we get a ~ frz , b ~ t~z and p = which is the well known 
vacuum solution of Chodos and Detweiler [21 J . 
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(iii) Further for c = a<i = 0, we recover the isotropic solution of Gr0n [22]. 
If we now calculate the 3D deceleration parameter, 

a z /a z 3 

such that k > 3 for accelerating model. From equation (12) it is evident that as we are 
discussing an expanding universe k < 6. Moreover, the second metric co-efficient b is 
regular when k < §. Further, when c > and < k < 3 we get dimensional reduction 
but no acceleration. It also follows from equation (13) that for (ai — cr 2 ) > 0, c > and 
< k < 3, the mass density also remains non- negative as is also evident from figure 1. 
Let us consider the situation for c < 0. In this case p is always positive. If we further 
subdivide the situation for (i) < k < 3 and (ii) 3 < k < |, we see for case (i) p > 0, but 
no dimensional reduction is possible and no acceleration, for the case (ii) while p is again 
positive, we interestingly get both dimensional reduction and acceleration. However small 
extra dimensions have of late been somewhat out of favour following the resurgence of 
interests in different brane inspired cosmological models where the extra dimensions need 
not be small to account for the vexed hierarchial problem of particle field theory. 

It is conjectured that during the process of dimensional reduction the extra dimensions 
finally stabilize at a very small length and loose their dynamical character before the extra 
scale factor b vanishes. Thereafter the cosmology enters the 4D phase without having any 
reference to the extra dimensions [29] . So the second singularity at b = is never reached 
and the invariants of equations (16-17) as also the mass density never have the chance 
to diverge there. For this model this transition has far reaching implications because 
the very existence of the extra space, so to speak, seems to induce inhomogeneity in 
this case. So not only do we enter a 4D era, it also envisages a smooth transition from 
a multidimensional, inhomogeneous phase to a 4D homogeneous one. Interestingly this 
transition takes place without forcing ourselves to choose very special initial conditions as 
is the practice in the four dimensional cosmology. This, in our opinion is a very important 
aspect of our model. So the primordial inhomogeneities die down in a natural way. 



3. Homogeneous and Isotropic Universe 

As our universe is now manifestly homogeneous and isotropic we assume a = b and 
also the metric coefficients depend on time alone. These assumptions also facilitate to 
make comparison with the observational results as also to gauge the influence of extra 
dimensions on the recent observational findings. 

With these two assumptions our general field equations (1) now reduce to 
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Figure 1: The time evolution of a, b, p and A are shown in this figure for c > and < k < 3 
(say, (ai — cr 2 ) = 10, k = 2, a2 = and c = 1 In i/iis case dimensional reduction is possible 
but no acceleration. At compactification t = t c , p surges to infinity and cosmology enters J^D 
phase 



As there are two independent equations with three unknowns we once again assume 
A = kH 2 . With this assumption we finally get, 

2(rf+2) 
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we can set the constant C = 1 and -D = without loss of generality in equation (21). 
Further the positivity of p gives the restriction k < ( d+2 K li+3 ) 
As a special case to our solution (21) we see that 

2 

(i) k — ( i.e. , A = ) implies a = t d + 3 which is the higher dimensional generalization 
of the well known Einstein- de Sitter solution. 

(ii) further for the usual 4D case (d = 0) a ~ £3 which is the FRW matter dominated 
case. 

2 

(iii) on the other hand d = 0, but A^O , a~ t 3 - fc 
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Characteristics of the Model 



From equation (19) 

P + A _ i (05) 

or Q m +Qx = 1 ( in the absence of curvature term ) where Q m is the higher dimensional 
cosmic matter density parameter and f^A is the corresponding vacuum energy density 
parameter. From equations (21) - (25) it also follows that Vt m = ^+2)(d+3) 2fc anc ^ ^ A = 
(d+2)(d+3) suc ^ ^at ^ m + = 1- So our solutions are consistent. It also follows from the 
expression of ^a that the arbitrary constant k is also a measure of the cosmic vacuum 
density parameter. We further see that the deceleration parameter comes out to be 

a (d+l)(d+2) -2k (d+ 1) - (d + 3)Q A 
q = -^lP = 2{dT2) = 2 (26) 

For our cosmology to accelerate k > g±!K^±g) sucn that we finally restrict k as 
(rf+i)(rf+2) < ^ < (d+2)(d+3) ^ This j ne q Ua }ity ensures that we get both positivity of p 
and acceleration. To make our analysis consistent with the present day observational 
results we find that for q = —0.5, k = 2 • ^or ^ ne particular case of 5D universe 
(d — 1), taking the most acceptable value of H = 72 km/s/mpc we calculate the age 
of the universe to be equal to 25.5 Gyr, quite large. Almost same age was obtained by 
Viswakarma [51] in a different context. Another point to mention is that the above value 
of k implies (for d — 1) Q m o = 0.25 and ^ao = 0.75, which are fairly consistent with 
current observational results. 

Another attractive feature of our model is that it gives a sufficiently large age of the 
Universe as mentioned earlier. It is evident from our equation (24) that the age of the 
universe for our model comes out to be to = ^ e _^ Ho (for d — 1). This is very remarkable 
in view of the fact that the age of the universe in the FRW model with a constant A is 
uncomfortably close to the age of the globular clusters t GC = 12.5 ± 1.2 Gyr [25]. The 
quintessential models give even lower age. As there is considerable uncertainty in the 
exact value of H we have calculated ^ao and t for different value of k, keeping H a 
constant for a particular table as shown below. 

Next we have plotted the age of the universe to against f^Ao on the basis of our model. 
It is found that as f^Ao increases to also increases. If the required mass density Q m o was 
smaller or ^ao is larger one could get higher age in these model as is clear from figure 2. 

Again we see from the table 1 that the gradual increases of H provide the more 
realistic age of the Universe, for H = 100 km/s/mpc, the age of the universe to — 18 
Gyr whereas the most acceptable age of being 14 Gyr [53]. We also see from the table 1 
that for k = 3, q = 0, i.e. , the flip occurs for this value of k. For H = 100 km/s/mpc, 
the flip time £/ = 9 Gyr, i.e., before it the universe decelerates. So our model A ~ H 2 is 
amenable to past deceleration and current acceleration. 
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Table 1: Age of the universe for different values of Hq (for d = 1) 



k 


q 






t(Gyr) 


t(Gyr) 


t(Gyr) 










H = 72 km/s/mpc 


H = 80 km/s/mpc 


H — 100 km/s/mpc 


0.00 


1.00 


0.00 


1.00 


6.25 


5.63 


4.50 


1.00 


0.67 


0.17 


0.83 


7.50 


6.75 


5.40 


2.00 


0.33 


0.33 


0.67 


9.38 


8.44 


6.75 


3.00 


0.00 


0.50 


0.50 


12.50 


11.25 


9.00 


3.50 


-0.17 


6.75 


8.44 


11.25 


13.50 


16.88 


4.00 


-0.33 


0.67 


0.33 


18.75 


16.88 


13.50 


4.50 


-0.50 


0.75 


0.25 


25.00 


22.50 


18.00 



4. Influence of Dimension 

(i) Age of the Universe : 

If we calculate t = h ( see e q ua ti° n (24)), taking the most acceptable values 

of f2 m0 = 0.33 and H = 72 ± 8 km/s/mpc and then plot t ~ d curve the age of the 
universe is found to decrease with the number of dimensions of the universe. This is a 
very interesting result of our whole analysis ( see figure 3). 

(ii) Cosmological constant : 

The relation between A with the number of dimensions of the universe is written 

as 

A = I# 2 (l-a m0 )(d + 2)(d + 3) (27) 

If we plot A ~ d one sees that A increases as the number of extra dimensions (d ) increases 
(figure 4). 

(iii) Deceleration Parameter : 

q = - [(fl m0 d - 3(1 - Q mQ )} (28) 

The above equation implies that as the number of dimensions increases q increases. 
For constant Q m o it is shown that q oc d. The figure 5 shows that (for Q m o = 0.33 ) as 
the number of dimensions increases, q increases. Obviously the flip is delayed with the 
increase of number of dimensions for any particular value of Q m o- 

5. Some Astrophysical Parameters in our model : 

In this section we shall very briefly discuss the proper distance D(z), luminosity 
distance (-Dl) and look back time. As these quantities are extensively discussed for 4D 
space time in the standard text books [26] we shall omit detailed mathematical calcula- 
tions and give only the final results generalised to ( 4 + d ) dimensions for our dust model. 
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Figure 2: As Q\q increases ( or Q m Q decreases ) to increases 



(i) Proper Distance : Let a star located at r = ro emit radiation at t = t\ and an 
observer at r = receive the same at t = to- The path being evidently a null geodesic 

(#1,02,03 , 0n — constants). The proper distance between the source and observer is 

given by 



D{z) = a 
for small z it reduces to, 



a °- = H.- 1 
ad 
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1 - (1 + Z) ^+2) 
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[l + q)z 2 + 



(30) 



Again, for 2 — > 00 but no acceleration i.e., k < \{d + i)(d + 2) 



D(2 = OO) = #n _1 



2(rf+2) 



(d+l)(d+2)-2fe 



2n 1 , for acceleration obviously -D(-z) — > 00 



(ii) Luminosity Distance : Luminosity distance is another very useful concept 
of relativistic astrophysics. If Dl is the luminosity distance of the object, L being the 
total energy emitted by a Galaxy in unit time then 



D, 



-) 

Aixl ) 



(31) 



After some straight forward calculation we get for our model 

2(rf + 2) 



H D L 



(d+l)(d + 2) -2k 



1 + z) 



2fc-(d+l)(d+2) 
1 + Z) 2(d+2l 



(32) 



10 



28 




Figure 3: This figure shows that as the number of dimension increases age of the universe 
decreases 



for small z, 



HnD, 



z + ^[5-Q m0 {d+l)]z 2 



(33) 



We have plotted Di{qo, 

Zj clS db function of z for various parametric values of go i n the 
figure 6. Note that all curves start off with the linear Hubble law z = H Q Di for small 
value of z, but then fan out, with only the curve for q = 1 staying linear all the way. As a 
rule we notice that, for the same redshift, the luminosity distance is larger for lower values 
of q . Thus for q = 1, we have, D L = -jr whereas for q = 0, we get D L = -^-[1 + 



again for q = —0.5, D L 



Ho 



(iii) Look back Time : The time in the past at which the light we now receive from 
a distant object was emitted is called the look back time. The radiation travel time ( or 
look back time ) (t — to) generalised to higher dimensions for photon emitted by a source 
at instant t and received at to is given by, 



2(d+2) 

In our case, a = t ( d+2 )( d + 3 )- 2 * so, 



H (t - 1) 



t-tr 



2(d + 2) 



ao da 
a 



(34) 



(d + 2)(d + 3) -2k L 
For small z the above equation reduces to, 



(d+2)(d+3)-2k 
1-(1+Z) 2 ( d + 2 ) 



(35) 



H Q (t - 1) 



(d + 2)(d + 3) -2k 2 
Z 2(d + 2^ Z + 



i - -(d + 3)a m z + 



(36) 
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Figure 4: As the number of dimension increases A also increases 



It follows from above that for a fixed z the look back time decreases with number of 
dimensions as evident from figure 7. Again it is shown that the look back time also in- 
creases with red shift (for small value of z). 



6. Cosmological model with constant A : 

For the sake of completeness as also of relevance to the current accelerating scenario 
we briefly discuss in this section an earlier work of one of us with a constant A [7]. In this 
work we discussed first zero pressure inhomogeneous 5D model with a constant A and 
later generalized it to a fluid obeying an equation of state p = mp ^ p 5 where p is the 
isotropic three pressure, p§ is the pressure in the fifth dimension, p is the matter density 
and m is a constant. At this stage it will not be out of place to point out that our set of 
solutions in Section 2 does not automatically reduce to the above mentioned works with 
a constant A as our ansatz A = kH 2 in the present work does not permit that type of 
reduction. 

For the line element (1) with d=l ( 5D space time) and A = constant we get for 
the matter dominated case (p = Pd = 0) 

a = ao sinh 5 pt (37) 
(b — cr 2 ) sinhpt + bi cosh.pt + % 



sinh 2 pt 



(38) 



p = J ^y-, where (do, bo, b\ are arbitrary constants, and c is the usual 4D curvature in the 
t-constant hyper surf ace. 
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Figure 5: q increases with the number of dimensions 



To facilitate comparison with our homogeneous space of section (3) we put c = such 
that the line element read 

i b sinh pt + bt cosh pt 
a = do smh 2 pt; b = j (39) 

sinha pt 

Following [27] we, at this stage, define a deceleration parameter such that an affective 
'scale factor' L becomes 

L = (a 3 b)* = sinh* pt(sinhpt + cosh pt)* (40) 

( putting ao, bo, bi unity, for simplicity). The well known 4D deceleration parameter ( say 
qi ) in this formalism becomes 

qi = — : — = —8 sinhpt cosh pt + 8 cosh 2 pt — 5 (41) 
L 2 

The behaviour of q\ here is interesting. At the early epoch the last expression shows that 
qi > ( good news for structure formations ). 

The flip occurs at the instant t = to, then the universe starts acceleration in conformity 
with the present day observation. It should, however, be emphasized that any attempt 
to compare the flip time of our model given above with observational findings will be 
too ambitions due to several initial assumptions in solving the equations as also due to 
the presence of so many arbitrary constants. One may, however, fine tune the arbitrary 
constants to account for the current observations. 

If we further assume b\ = we get isotropic expansion , i.e. , a = b. In this case the 
expression for q ( say 02) turns out to be 

cosh 2 pt — 2 
cosh 2 pt 
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Red Shift (z) 



Figure 6: The Luminosity distance Dl expressed as a function of the redshift z for q = 
1,0,-0.5. The relationship is linear, as predicted by Hubble's linear Law. For qo = —0.5, Dl 
increases faster. All curves merge for small z 



From the expression (41) and (42) it is evident that the flip occurs earlier for the anisotropic 
case, which is also evident from our figure 8. However the cosmological implications of 
this finding need further investigation. 

To end the section let us find out if the presence of extra dimensions, in any way, 
influences the instant of flip for the isotropic case. From the divergence equation in 
(3+d+l)-dimensional homogeneous space time for line element (1), we get, 

d b \( p + p ) = (13) 




b 

which for our case (p = 0, a = b) yields 

p = Ca- {d+3) (44) 
Using equation (44) in equation (19), we get 

Skipping mathematical details and adjusting the arbitrary constant C we finally get from 
equation (45 ) 

2 

a = a sinhd+3 pt (46) 



where P = y f+f-A-. So the deceleration parameter reduces to 

(d + 3) -2 cosh 2 pt 



(47) 
2 cosh pt 
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Figure 7: Look back time increases with red shift (for small value of z ) 



so the flip occurs at 



cosh pt c = 




Obviously the flip is delayed with the increasing of number of dimensions although a 
physical explanation of this behaviour is beyond the scope of the present investigations. 
Interestingly we got the same result in section 4 for the case of varying A ( figure 9). 

7. Discussion 

In this work we have studied a multidimensional model (both homogeneous and inho- 
mogeneous) with a time varying cosmological parameter. It is found that under certain 
conditions the cosmology gives an accelerated expansion. In the inhomogeneous case no 
initial deceleration is found. However for the homogeneous case an initially decelerating 
universe starts accelerating undergoing a flip. This is good for both structure formation 
and current observational status. Some of our cases are also amenable to the desirable 
feature of dimensional reduction. As mentioned earlier the dimensional reduction for 
our inhomogeneous model seemingly points to the fact that a primordial inhomogeneous 
higher dimensional cosmology enters a 4D homogeneous era although we do not have 
to make any stringent initial conditions to achieve this homogenization as tn the stan- 
dard 4D case. This is an important finding of our analysis. While it is conjectured that 
the extra space stabilises at planckian length it is not apparent from our analysis how 
that mechanism in the form of a sort of repulsive potential field actually works in our 
model. This is definitely a defect of our model. However in an earlier work Guendelmann 
and Kaganovich[28j studied the Wheeler- Dewitt equation in the presence of a negative 
cosmological constant and dust and showed that quantum effects do stabilise the volume 
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Figure 8: The time evolution (t not to scale) of q is shown in this figure. It is shown that flip 
for qi is earlier than flip for q<i . 



of the universe, thus providing a mechanism of quantum avoidance of the singularity. 
Other two important findings in our analysis are that the time of flip depends on the 
number of extra dimensions. In fact as the number of dimensions increase the flip is de- 
layed thus higher dimensional models allow more time for structure formation. Secondly 
the flip also depends on anisotropy. Shear simply hastens the arrival of the instant of flip. 
It is too premature to attempt a physical explanation of the above two results but the 
results are interesting enough to warrant further investigations in this regard. To end a 
final remark is in order. We mentioned in the introduction that both extra dimensions 
and inhomogeneities produce a set of extra terms to the 4D FRW model, which creates a 
sort of back reaction [12] in the process. While the effect of extra dimension is discussed 
here in some detail, the influence of inhomogeneities in higher dimensional model will be 
discussed in our future work. 

To end a final remark may be in order. We have here presented results which include 
cases when the extra dimensions also inflate along with the ordinary ones in section 3. 
While the idea of large extra dimensions is not as repugnant these days as in the past 
following its new found relevance in brane models but an inflating extra dimension goes 
far beyond that. This is a serious shortcoming of the present analysis, particularly the form 
of the spacetime we chose in section 3. While working in higher dimensional spacetime 
one should see that the extra space should form a compact manifold with the symmetry 
group G so that the (d + 3) spatial symmetry group is a direct product 0(3) xG and not 
0(d + 3), although there exists a good number of works in this type of space also. 

Acknowledgment : We acknowledge the financial support of UGC, New Delhi as 
also valuable comments of the anonymous referee. 
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t (not to scale )- 



Figure 9: Flip is delayed with the number of dimensions 
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